Random walk weakly attracted to a wall 
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Abstract: We consider a random walk X n in Z+, starting at Xq = x > 0, with transition 
probabilities 

P(X n+1 = X n ± l\X n = y > 1) = - =f 

V n+l n in y _ ; 2 4y + 25 

and X n+ i = 1 whenever X n = 0. We prove EX n ~ const, n 1- ^ as n / oo when 5 G (1, 2). 
The proof is based upon the Karlin-McGregor spectral representation, which is made 
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CN ' explicit for this random walk. 

KEYWORDS: Random walk, orthogonal polynomials, pinning, wetting 
AMS subject classification: 60J10, 82B41, 42C05 



C3 



1. Introduction 

Random walks have been used in many different fields of physics, economics, biology... 
Usually, it evolves in a translation invariant environment or a random environment whose 
average is translation invariant. Random walks in inhomogeneous environment, in particu- 
lar with a reflecting or attracting or repelling wall, have been used to mimic the behaviour 
of a liquid interface on top of a solid substrate, wetting phenomena, and a variety of other 
phenomena associated with the pinning of an interface. This is the motivation for the 
model presented below, whose specific form was chosen so as to derive rigorously and ex- 
plicitly the result of competition between a relatively long ranged attraction and reflection 



(N 



O ' at the wall. 

We consider a random walk X n in Z + , defined by X = x > 0, 

OO ' 



O. Py =F(X n+1 =X n + l\X n = y > 1) 



1 6 y 



2 4y + 2d 2y + d 

><: «y = ^ x ^ = Xn-i\x n = y >i) = \ + ^ = ^ 5 

and X n+ i = 1 whenever X n = 0, i.e. po = 1, qo = 0. Let P x (-J denote the corresponding 
probabilities, and E x (-) the corresponding expectation values. 

The walk obeys the detailed balance condition with respect to the measure 7r on Z + 
defined up to a multiplicative constant by 

w tt p. „ (2y + <s)r(<s + i)r(;,) 
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which obeys n y ~ y~ when y — > oo. We restrict our attention to 5 > 1, and normalise 7r 
as a probability measure, with 7ro = (5 — l)/(2<5). The dynamics is periodic, with the odd 
and even components exchanged under one step of the dynamics, tv = (ir even + 7r dd)/2, 
with /Teven = (27To, 0, 2n 2 , . . .), and 7r dd = (0, 2tt\, 0, 27T3, . . .). Starting from Xq = 0, we 
have convergence in law X 2n — > X^ en and X 2n +i — > X^ d . The first moment exists only 
for 6 > 2, with EX X = E(A^ en + X£ d )/2 = 6/(2(6 - 2)). 

We focus our attention to 1 < 6 < 2, so that E x X n — > oo as n / oo. For |z| < 1 let 
g e (z) = ^oX n , 9o (z) = ^ z n E X n (1.3) 

n even n odd 

Our main result is the following theorem: 

Theorem 1: Let 6 G (1, 2). Then as z / 1 or n /" oo, 

= r(2 - (i - ,)-" 2 (i + - z)^( s -^ (i.4) 

E X n = Ksn 1 -? (l + o(n-< 1 -M'-»y) (1.5) 
wiiere is a constant depending upon 6 as defined in (4.12). 

Remark 1: The constant K$ /* +00 when 5 \ 1 or 6 S 2, but the O(-) are not uniform 
in 6, the theorem says nothing about 6 = 1 or 6 = 2. 

Remark 2: The odd generating function g (z) obeys the same bounds as g e (z). For any 
starting point x, E x X n obeys the same bounds as E X n . 

Remark 3: (1.5) follows from (1.4) by a Tauberian theorem [F, Thm. 5 p 447]. This uses 
monotonicity of E X n for n even, which holds for any reflected random walk, as recalled 
below together with more specific monotonicity arguments. 

Consider two walks: X n started at x, with parameter 6, and X' n started at x', with 
parameter 6' . Assume 6 > 6' and x < x' and x' — x even. The two walks may be coupled, 
e.g. using a single random number for both walks when they meet. This implies that 
E x X n is for each n an increasing function of x over the even integers, and also over the 
odd integers, and a decreasing function of 6 G (—1, +00). Monotonicity in 6 is an example 
of a general monotonicity property in the transition probabilities p y , q y . Then 

F (X 2 = 0)E X n + F (X 2 = 2)E 2 X n > E X n (1.6) 

Therefore E X n is an increasing function of n over the even integers, and also over the odd 
integers. And the same property holds for EiX n . But E x X n for x > 2 is not monotonous 
in n, even or odd. We have 

E(X n+1 -X n \X n = y) = --^— , y>\ (1.7) 

2y + 6 
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so that E X 2 n < E X 2n -i- And since E X n is increasing over the even (or the odd) 
integers, we have E X 2n +i > E X 2n . 

(1.4) will be proven in Section 4. In Section 2 we exhibit the orthogonality measure of 
our random walk polynomials: this is Theorem 2, our main technical result, opening the 
way for Theorem 1. In Section 3 we solve a differential equation for generating functions, 
differential equation associated with the recursion formula for random walk polynomials. 

Our random walk is of a very special form, but it could be used for comparison or 
as input for more realistic models. For example the polymer pinning model of Alexander 
and Zygouras [AZ] combines i.i.d. disorder with a spatially inhomogeneous Markov chain, 
which could be built from our random walk. 

Other aspects of this model and closely related models will be treated in a forthcoming 
paper [H]. 

2. Random walk polynomials and their orthogonality measure 

Our first tool is the Karlin-McGregor representation theorem [KM]: let L 2 {ir) denote 
the Hilbert space of complex sequences (f y ) ye z + obeying Xl^Lo \fy\ 2rK y < 00 • Then 

(Tf) y =p y f y+1 + q y f y - 1 (2.1) 

defines in L 2 (-k) a self-adjoint operator T of norm less or equal to one. Let eo = (1, 0, 0, . . .). 
The Karlin-McGregor representation theorem gives 

F x (X n = y) = ^(T n Q x (T)Q y (T)eo,e ) (2.2) 

where {Q y (t)} ye j j+ is a family of polynomials of degree y in t, defined recursively by 
Q-i = 0, Qo = 1, and 

tQy = PyQy+l + q y Qy-l , V > (2.3) 

giving polynomials Q y of degree y and parity (— l) y , with Q y (l) = 1 for all y > 0. Using 
the spectral resolution {E t } of the self-adjoint operator T, and d/j,(t) = d(E t e 0l e ) : one 
gets 

F x (X n = y)= 7 ^ [ d/j(t)t n Q x (t)Q y (t) (2.4) 

which implies that {Q y (t)} ye z + is a family of orthogonal polynomials in the probability 
measure dfj,(t), 

[ dn(t)Q x (t)Q y (t) = ^5 x , y , x,y>0 (2.5) 

7-1 TTy 

Given the family {Q y {t)} ye j j+ , (2.5) characterizes a unique probability measure dfi, termed 
the orthogonality measure of the family. Letting n / oo in (2.4) shows that [KM, pp 70-71] 

dn(t) = ir (s(t - 1) + S(t + lfj + d/j c (t) (2.6) 
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where dfi c (t) is absolutely continuous with respect to the Lebesgue measure. Indeed as 
n y oo with n + y — x even, the LHS of (2.4) tends to 2ir y , while the contribution from 
dfi c to the RHS tends to zero. 

Our second tool is Dette's theorem [D]: the orthogonality measure dfi 1 of the first 
associated polynomials of a random walk is related to the orthogonality measure dfi* of 
the dual random walk through 

d fJ , 1 (t) = —{l-t 2 )dfj,*{t) (2.7) 

The first associated polynomials Q y are defined by Qli = 0, Q\ = 1 and 

tQ^Py+iQl+i + qy+iQl-i, y>o (2.8) 

The dual random walk polynomials Q* are defined by p* — q y and q* = p y for all y except 
Pq = 1 and = 0, so that Qq = 1, Q\ = t and 

tQ^piQi+i + qlQl-^qyQl+i+PvQl-!, V>1 (2-9) 
Dette's theorem may be applied starting from the dual, giving 

d/JL*> 1 (t) = —(l-t 2 )diJL(t) (2.10) 
Pi 

The first associated dual polynomials Q*' 1 are defined by Q*l\ = 0, Q^ 1 = 1 and 

tQ*/ = ?tf+ iQSi + » ^ ( 2 - n ) 

The general definitions (2.3) (2.8) (2.9) (2.11) are made explicit by inserting our p y } s and 
Qv% giving 

(2y + S)tQ y = yQ y+1 + (y + d)Q y - 1 , y>l, Qi(t)=t 
(2y + 2 + S)tQl = (y + + (y + 1 + S)Q 1 y _ 1 , y>0 

(2y + s)tQ; = ( y + s)Q; +1 + yQ;_ 11 y >i (2 ' 12) 

(2y + 2 + 5)tQ*/ = (y+l + 5)Q^ 1 + (y + l)Q* y \ , y > 

We thus have four distinct families of orthogonal polynomials, and the corresponding four 
distinct orthogonality measures fi, fj}, fx*, fx*' 1 . 

The third and last step is to relate our polynomials to Gegenbauer polynomials of 
index A, defined by = 0, Gq = 1 and 

(2y + 2A)tGj = (y+l)G} +1 + (y-l + 2A)Gj_ 1 , y > (2.13) 
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It appears that Q y = G§ +1 , so that Q*'* = Gy +1 '* , and also Q* = p^|^ . These in 
turn satisfy 

| +M _ (y + l)ir(*+l) 1,1 _ (y+l)£(g5) i^.y 

r( y + i + 5) Gy r( y + ^) Py M) (2 - 14) 



where P"'^(t;c) are c-associated Jacobi polynomials, or Wimp polynomials [W]. The or- 
thogonality measure da** 1 is therefore aJ 
polynomials, namely [W, Th. 3 p. 996] 



S-l S-l 

thogonality measure da*' 1 is therefore also the orthogonality measure of the P y 2 ' 2 (£; 1) 
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d/i*' 1 ^) = — — dt I normalisation (2-15) 
\F(t)\ I 



where 



m - ,* (u - ^ + (4^ ^ ¥^ T) 

= -l + 0((l-t)^) as 
with 2-F1 the Gauss hypergeo metric function and 



(2.16) 



r(*)r(±=*) 



if — v V~r (2.17) 

r(*=i) 



We have F(-£) = -F(£)*. Then (2.6)(2.10) and (2.15) give 



(l-£ 2 )^ , 

d/j, c (t) = — — (it / normalisation (2.18) 
\F(t)\ I 



which yields: 



Theorem 2: The orthogonality measure of the family of polynomials defined by (2.3) 
with (1.1) and 1 < S < 2 is the even probability measure on [—1, 1] defined by 

dfji(t) = iro(6(t - 1) + 6(t + if) + dfi c (t) (2.19) 

where S(-) is the Dirac measure at and d\i c is given by (2-16)-(2.18) and tt = with 
the normalisation f_ dji c = 5~ l . 
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3. Generating function 

Using (2.4) and X n < X + n, we have 



n „i 

E X n = V^/ dii{t)t n Q y {t) (3.1) 

For n even, using Q y orthogonal to Qo = 1, and for y odd Q y also orthogonal to t n , and 
then using (2.6), we have 



n „i 
E X n = V / cfc c (f)(f n - l)Q y (f) 

even 

«1 n 

= 2 / dp c (f)(f»-l)£^Q tf (f) 

Jo ^0 
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The corresponding generating function is defined as 

g e (z) = z n E X n 

n even 

neven v=o 

even 

u y even 

For n odd, using Qj, orthogonal to Qi = t for y > 2, and (2.5) for y = 1, 

n „i „i 

y~{ 7T0 y_l 7T0 J-l 

= 1 + 2 / d// c (t)(t n - t) V ^Q„(t) = E 1 X n . 1 
Jo U tto 

odd 

and the corresponding generating function 

9o (z) = z n ^oX n 

n odd 

= -r^— 2 + 2 f df(t) ~ t) E 

^ nodd ,=i ^ 

odd 

y odd 
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(3.2) 



(3.3) 



(3.4) 



(3.5) 



The recursion (2.3) (2.12) defining the random walk polynomials gives a poor uniform 
bound for these polynomials, e.g. \Q y (t)\ < 3 y for t e (—1, 1). We define, for \u\ < 1/3, 



oo 



= E ^U(*K = r(5 + 1) f> + 8) T( T +f +1) Qv(t)u v (3.6) 
and aim at an analytic continuation giving, for |z| < 1, 

(3.7) 



1 - z 2 t 2 l-z 2 J V 1 - 1 - z 2 

%{zt) %{z)\ , f%(-zt) %(-z) 



1 - z 2 t 2 1 - z 2 

The recursion (2.3) (2.12) may be converted into a first order differential equation for ^ (w). 
We first get a differential equation for 

oo -p/ \ oo 

*«(«) = T(S + r(v + y 5 + l) Q y {t)uy = E H v^ uV ( 3 - 8 ) 

y=l ^ + + ^ y=l 

with 

obeying, with any arbitrary value for if , 

(2y + 8)tH y = (y + S + l)H y+1 + (y - l)H y _ x , y > 2 (3.10) 

with 

Hl{t) ~ ITT ' m) ~ (s + i)(s + 2) ~ (sTTW+2) (3 - n) 

Multiplying (3.10) by u y and summing over y > 2 yields 

(1 - 2tu + u 2 )& t (u) = t-u- 5(w _1 - t)$ t (u) (3.12) 
with $t(0) = 0, whose solution is 

$ t ( u ) = 8~ 1 - U - 5 (l-2tu + u 2 )i dvv 5 - 1 (l-2tv + v 2 )-^ (3.13) 

We thus get: 

Lemma 3: The function ip t (u) defined in (3.6) may be expressed as 

*,(«) = 2«^(«) + M>«(«) = -1 + j4^^2 ( 5 ~ 1 - ( 3 - 14 ) 



where 3>t( u ) JS the solution of the differential equation (3.12). It extends to an analytic 
function in the disc \u\ < 1. Its derivative may be expressed as 

w;M= /c-" 2 ) mi r , (3 . 16 ) 

u(l-2tu + u 2 ) u s (l-2tu + u 2 ) 2 -Uo (l-2tv + v 2 )i 



with 



B t (u) = 4u(l - £) - 2t(l - w) 2 + ^ Ut){1 ^ (3.16) 



u 



4. Proof of (1.4) 

The leading order and next to leading order in (3.7) as z / 1 are found in 

f 1 , . , fzt^'Azt) z%(z)\ /A ^ 

Indeed the leading orders come from the singularity at £ = 1 in the integral. In (3.7) 
there is a symmetry or anti-symmetry as z — > — z and £ — > — £ jointly, associated with 
the even/odd symmetries. The terms not included in (4.1) correspond to — z ~ —1, not 
singular with £ > 0. For 1 - z C 1 and 1 - i <C 1 we have 



£ 1 „ l + 2 2 £ 2 ll-£ 



-(l-£) 



1-zH 2 1-z 2 v '(l-z 2 * 2 )(l-z 2 ) 21- z(l- z) + (l-t) 
t 1 _ £(1-^ 2 ) + l-z 2 t 2 1 1 2(l-z) + l-£ 



(4.2) 



1-zH 2 1-z 2 (1 - z 2 t 2 )(l - z 2 ) 21-z(l-z) + (l-t) 
(3.15) (3.16) may be written as 

qt { z )H' t { z ) = 5^—^q t {z)-B t {z)5z- 5 q t {z)i f dw^q^v)^ (4.3) 
z Jo 

with 

q t (z) = l-2tz + z 2 = (l-z) 2 + 2z(l-t), q' t (z) = 2(z - t) (4.4) 
Together with (4.2), we have to estimate %(zt) ± %(z) as z / 1, £ /" 1. We have 

l-z 2 £ 2 1-z 2 /. l-£ 



+ («) 



z£ * 

B t (z) ~ 4(1 - £) - 2(2 - 8)(1 -z)(l-t) + 2(5 - 1)(1 - z) 2 

~ 2(2 - 6){1 - t) - 4(5 - 1)(1 - z) (4 6) 

S t (zt)~S t (z)(l-(l-t)^) 



r 1 dw 5 - 1 _ r 1 

Jo q t {v)i Jo ((i- v 



dv 



((l-t/)' + 2(l-f))* 
dx 



+ / dv 



,5-1 



(l-f)T- / 

JO 



(x 2 + 2)' 



Qt(v)i ((i_ v )2 + 2 (i_t)) 
+ 0{1) 



s _ 2-^0Fr(^i) 
r(I) 



= (i-t)- 



The range (0, 1) of the integral in (4.1) is split into four intervals, according to 
< 1 - (1 - z) a < 1 - (1 - zf < 1 - (1 - z) 1 < 1 



(4.7) 



(4.8) 



with 0<a<l</3<2<7 and respective contributions denoted gf , gf 13 , gf 7 , gj. We 
begin with the leading contribution, gf . 

• I- (I- z) a <t <l- (I- zf : 
q t (z) ~ 2(1 - t) and B t (z) ~ 4(1 - t). Then (4.3) with (4.7) yields 



(1 - tf%{z) ~ 5(1 - Z )(l - t) - (1 - f)*^ 



r(|) 



(4.9) 



and 



l-(l-zf , . 

rft(l-t 2 ) — 



-(!-*)< 



1 - Z 2 t 2 1 - Z S 
S-l \ f l-(l-z) f: 



2 r(f) 



2) Ji-d-zy 



s-s pK r( 



2—5 



(5-1' 



2 r(§) 



■(l-z)i 



(i- 2 ) c 



■ -2 1 _ * 



(l_ z )/3-l 1 + X 



l-z(l-z) + (l-t) 



= 2 ^ 5 \/f r (~y~) r ( 1 - 1) a - -) | - 2 ( 1 + - *) (1 -° )(1 -*> + (i - 

(4.10) 

Collecting all previous error terms we get 



a/3 

9i 



with 



( z ) = r(2 - - z)i~ 2 [1 + - z)* 1 "")* 1 -*) + (i _ *)«^ + (1 - ^) J 



= 2 — 



2 1-1 / y_! \F(t)\ 2 



-dt 



(4.11) 
(4.12) 



and F(t) given by (2.16). The five error terms in 0(- ■ •) in (4.11) come respectively from: 

1) Extending the x-integral to +oo in (4.10) 

2) F{t) = -l + 0((l-t) !L ?) 

3) ^f' t (zt) — ^t(z) ~ 5, derived from (4.3)-(4.6), with the leading contribution from (4.5). 

4) Extending the x-integral to in (4.10) 

5) q t (z) ~ 2(1 - t) (l + ^) and B t (z) ~ 4(1 -t)(l + (S- 1) |=^g) . 
Choosing 

a= sup fniin{(l-a)(l--),a— ^— }) =2-5 (4.13) 
5e(i,2)^ 2 2 J/ 

and /3 within 4 - 2^ < (3 < 2, e.g. /3 = 3/2, then brings (4.11) to (1.4). The remaining 
ranges in the t-integral also fall within the previous error estimates: 

• < t < 1 - (1 - z) a : 

g « = (i - z )i- 2 o((i - ^a-^a-D + (i _ z )t + (i _ 2)-^) (4.i4) 

• 1 - (1 -zf < t < 1 - (1 -z)"<: 

gp = (1 - z)!- 2 o{il - z)W-»i + (1 - z) 2 ^) (4.15) 

• 1 - (1 -z)i < t < 1: 

^ = 0((l-z)- 1 ) (4.16) 
This completes the proof of (1.4) and Theorem 1. 
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